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Abstract 

We provide new examples of diffusion operators in dimension 2 and 3 which 
have orthogonal polynomials as eigenvectors. Their construction rely on the finite 
subgroups of 0(3) and their invariant polynomials. 
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1 Introduction 

We investigate in this paper new examples of bounded domains in on which there 
exists a probability measure fi with an orthonormal basis of such that the el¬ 

ements of this basis are eigenvectors of a diffusion operator. To determine such a 
basis, one needs first to define a valuation (a definition for the degree) for a polyno¬ 
mial in two variables. The complete determination of all the possible such domains 
in has been carried in [3], under the restriction that the valuation is the usual 
one (that is the degree of the monomial is p + q). We shall show in this paper 
that relaxing this requirement on the valuation leads to many new models. We have 
no claim to exhaustivity, and for the moment have no clue about a possible scheme 
which would lead to a complete classification for the general valuation. However, 
the domains that we describe here all share some common algebraic properties that 
we want to underline. 

The construction of these domains rely mainly on the study of finite subgroups of 
0(3), and are in particular related to the Platonic polyhedra. It relies on the study 
of invariant polynomials for subgroups of 0(3). The analysis of these invariants also 
lead to the construction of new polynomial models in dimension 3. 
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2 Orthogonal polynomials and diffusion opera¬ 
tors 


The short description of diffusion operators that we present below is inspired from [T], 
and we refer the reader to it for further details. 

Diffusion operators are second order differential operators with no zero order 
terms, and are central in the study of diffusion processes, solutions of stochastic 
differential equations, Riemannian geometry, classical partial differential equations, 
potential theory, and many other areas. When they have smooth coefficients, they 
may be described in some open subset Q of as 

(2.1) L(/) = + J2b\x)di{f), 

ij i 


where the symmetric matrix is everywhere non negative (the operator L 

is said to be semi-elliptic). We are mainly interested here in the case where this 
operator is symmetric with respect to some probability measure /i, that is when, for 
any smooth functions /, ^, compactly supported in D, one has 

(2.2) [ fL{g)dii= f g'LOdii. 

Jq Jq 

We then say that /i is a reversible measure for L, which reffects the fact that, in a 
probabilistic context, the associated stochastic process has a law which is invariant 
under time reversal, provided that the law at time 0 of the process is /i. 

When fi has a smooth positive density p with respect to the Lebesgue measure, 
this symmetry property translates immediately in 

(2.3) h\x) = djg'-^ (x) + g^^dj log p, 

3 3 


which shows a fundamental relation between the coefficients of L and the measure /i, 
and allows in general to completely determine p up to some normalizing constant. 

Let us introduce the carre du champ operator T. For this, we suppose that we 
have in jC‘^{p) some dense algebra A of functions which is stable under the operator 
L, and contains the constant functions. Then, for (/, ^) E v 4 , we define 


(2.4) 


T{fJ) = -{L{fg)-fL{g)-gL{f)). 


If L is given by equation ( |2.1[ ), and when the elements of A are at least it 
turns out that 

r(/,5) = '^g'^difdjg, 




so that r describes in fact the second order part of L. The semi-ellipticity of L 
translates into the fact that r(/, /) > 0, for any f ^ A. If we apply formula (2.2) 
with ^ = 1, we observe that Lfdp = 0 for any f ^ A. Then, applying (2.2) again, 
we see immediately that, for any (/, ^) E A 


(2.5) 


[ fH9)dp = - [ T{f,g)dp, 
Jq Jq 
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so that the knowledge of F and p describes entirely the operator L. We call such a 
triple (f], r, /i) a Markov triple, although we should also add the algebra A. Thanks 
to ( |2.lD , we see that L(x^) = 6^ and r(x% x^) — . The operator F is called the co¬ 

metric, and in our system of coordinates is described by a matrix F = (F(x^,X'^)) = 

In our setting, we shall always chose A to be the set of polynomials. Under the 
conditions that we shall describe below, we may as well extend A to be the set of 
the restrictions to of the smooth functions defined in a neighborhood of fi, but 
this extension is useless in what follows. 

The fact that L is a second order differential operator translates into the change of 
variable formulas. Whenever / = (/i, • • • , /n) ^ A^, and whenever $(/i, • • • , /^) G 
M, for some smooth function $ : M, then 


(2.6) L($(/)) = 5] dMf)Ufi) + Yi fj 

i 

and also 


(2.7) 






When A is the algebra of polynomials, this applies in particular for any polynomial 
function $. Indeed, in this context, properties (2.6) and (2.7) are equivalent. 


As long as polynomials are concerned, it may be convenient to use complex 
coordinates. That is, for a pair of variables (x,^/), consider z = x + iy and z = 
X — iy^ using linearity and bilinearity to extend L and F to z and z, for example 
setting L(z) = L{x) + iL(y), V{z^ z) = F(x, x) — F(^, y) + 2iF(x, y). Then one may 
compute \j(yP{z^z)) and V[P{z^z)^Q{z^z)^ for any pair of polynomials P and Q 
in the variables {z^P) using the change of variable formulas ( |2.6[ ) and ( |2.7[ ). One 
may then come back to the original variables x and y setting x — {z + z)/2^ y — 
{z-z)/{2i). 

Moreover, we shall restrict our attention to the case where the matrix is 
everywhere positive definite, that is when the operator L is elliptic. In this situation, 
one may expect L to have a self adjoint extension (not unique in general), and then 
look for a spectral decomposition for this self adjoint extension. We may expect 
then that the spectrum is discrete, and look for the eigenvectors. 

It is quite rare that one may exhibit explicitly any eigenvalue or eigenvector, and 
this makes the analysis of such operators quite hard. However, a good situation is 
when there is a complete basis formed of polynomial eigenvectors, in which 

case one may have explicit computation for the eigenvalues and expect a good de¬ 
scription of the eigenvectors (recurrence formulas, generating functions, etc). These 
polynomials being eigenvectors of a symmetric operator are orthogonal whenever 
the eigenvalues are different, and this leads to a family of orthogonal polynomials 
for the invariant measure g. 

Unfortunately, this situation does not appear quite often. In dimension 1 for 
example, up to affine transformations, there are only 3 cases, corresponding to the 
Jacobi, Laguerre and Hermite polynomials, see for example [2]. 
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1 . 


The Hermite case corresponds to the Gaussian measure 
the Ornstein-Uhlenbeck operator 




dx on M and to 


^ d 

^OU = 3 -^ - ^ 3 -* 
dx^ dx 


The Hermite polynomial Hn of degree n satisfy YjojjPn = —nPn- 

2. The Laguerre polynomials operator correspond to the measure fia{dx) = CaX^~^e~^ dx 
on (0, 00 ) , a > 0, and to the Laguerre operator 


d? 

La = + (a - x) 


d 

dx 


The Laguerre polynomial with degree n satisfies = —nL^\ 

3. The Jacobi polynomials correspond to the measure fia^hi^x) — Ca^h{^~xY~^{1+ 
x)^~^ dx on (—1,1), a, 6 > 0 and to the Jacobi operator 


^2 

La,6 = (1 - - (a - 6 + (a + b)x) 

The Jacobi polynomial with degree n satisfy 

La,64“’''^ = -n(n + a + 6-l)jy''). 


d 

dx 


In this paper, we concentrate on probability measures on bounded domains C 
For such measures, the set of polynomials is dense in and we want to 

construct bases of formed with polynomials. There is not an unique choice 

for such a basis. 

First, we choose a valuation. That is, choosing some positive integers ai, • • • , a^, 
we decide that the degree of a monomial x^ • • • x^^^ is aipi + • • • a^Pd- Then, the 
degree of a polynomial is the maximum of the degrees of its monomials. 

This being done, for n G N, we look at the finite dimensional vector space Tin of 
polynomials with degrees less than n. One has T-Ln C T-Ln+h UnT-Ln is the vector 
space of polynomials. It is dense in Then, a polynomial basis is a choice, for 

any n, of an orthonormal basis is the orthogonal complement of Tin+i in 

Our problem is then to describe for which open bounded subsets C M^, one 
may find a probability measure p on it with positive density p{x) with respect to 
the Lebesgue measure, and an elliptic diffusion operator L on such that such a 
polynomial basis for p is made of eigenvectors for L, for some given valuation. We 
restrict our attention to those sets with piecewise smooth boundary. Let us call 
such a set a polynomial set, and the triple (i2,r,/i) a polynomial model. 

We recall here some of the results of [3], where the same structure is described 
only for the usual valuation (that is when all the integers ai are equal to 1), but 
easily extended to the general valuation case. We have 

Proposition 2.1. Choose a valuation deg described as above by some integer pa¬ 
rameters (ai, • • • ,a^); and let (i2,r,/i) be a polynomial model in M^. Then, with L 
described by equation 
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1. for i = l^ -'' b'^ is a polynomial with deg(6^) < a^; 

2. for i, j = 1, • • • ^d, is a polynomial with deg{g^^) < + aj; 

3. the boundary is included in the algebraic set {det{g^^) — 0}; 


4 . if Pi’’’ Pk = 0 is the reduced equation of the boundary (see remark 2.2 
below), then, for each q = 1, • • • k, each i = I, • - d, one has 


( 2 . 8 ) 


T(\og Pq, Xi) — 


where Li^q is a polynomial with deg{Li^q) < a^; 

5. all the measures with densities ^ak\Pi\^^ '' * \Pk\^^ on Q., where 

the OLi are such that the density is is integrable on are such that (f], F, 

is a polynomial model; 

6. when the degree of Pi - •• Pk is equal to the degree of det{g'^^) there are no other 
measures. 


Conversely, assume that some bounded domain Q is such that the boundary dQ 
is included in an algebraic surface and has reduced equation Pi - - - Pk =0. Assume 
moreover that there exists a matrix {g^^{x)) which is positive definite in Q and such 
that each component g'^^ {x) is a polynomial with degree at most aid-aj. LetT denote 
the associated carre du champ operator. Assume moreover that equation (2.8) is 
satisfied for any i = 1, - - - ,d and any q = 1, • • • ,k, with Li^q a polynomial with 
degree at most ai. 

Let {ai, - - - ,ak) be such that the \Pi\^^ - - - \Pk\^^ is integrable on with respect 
to the Lebesgue measure, and denote = Cai,--- ^a^Pi^ - --P^^dx, where 

Cai,---,ak normalizing constant such that is a probablity measure. 

Then (12, F ,is a polynomial model. 

Before giving a sketch of the proof of Proposition |2.H let us make a few remarks. 


Remark 2.2. We say that Pi - - - Pk = 0 is the reduced equation of dQ when 

1. The polynomials Pi are not proportional to each other. 

2. For i — 1, - - - k. Pi is an irreducible polynomial, both in the real and the complex 
field. 

3. For each i — 1, • • • ,k, there exists at least one regular point of the boundary 
512 such that Pi{x) = 0. 

4 . For each regular point x G 512, there exist a neighborhood V and of x and some 
i such that 512 n V = {Pi{x) = 0} H V. 

In particular, this does not mean that any point satisfying Pi{x) = 0 for some i 
belongs to 512. 


Remark 2.3. The determination of the polynomial models therefore amounts to the 
determination of the domains 12 with an algebraic boundary, with the property that 
the reduced equation of dLl is such that the set of equations (2.8) has a non trivial 
(and even positive definite concerning {g'^^)) solution, for g'^^ and Li^q. Looking at 
the form of these equations, given the reduced equation of dLt, they appear as a linear 
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homogeneous equation in the eoeffieients of the polynomials and of the polynomials 
Li^k- Unfortunately, there are in general mueh more equations to be satisfied that 
unknowns, and this requires very strong eonstraints on the polynomials appearing in 
the redueed equation of the boundary. 


Remark 2.4. The set of equations {2.8), whieh are eentral in the study of polynomial 
models, may be redueed to less equations, when k > 1. Indeed, if we set P — Pi • — P^, 
it reduees to 


(2.9) 


r(a:i,logP) = Li, deg(Lj) < ai. 


To see this, assume that this last equation holds with some polynomial Li. Then on 
the regular part of the boundary deseribed by {Pq{x) — 0}; we have T{xi,Pq) = 0^ 
sinee 


T{xi,Pq) = Pq{Li 


l^q 




Therefore, Pq divides V{xi,Pq) 


Proof. —(Of Proposition [24] ) . 

We shall be a bit sketchy in the details, all the arguments being borrowed 
from [3]. Let TLn be the finite dimensional vector space of polynomials P such 
that deg(P) < n. From the definition of a polynomial model, L{TLri) C TLri’ In 
the representation S of L, we have U = L(x^) and gU = T{xi,Xj). Therefore, 
U G TLai and, from the representation ( |2.4| ) of F, gU g TLai^aj- This gives items 
and [21 

Now, since L has polynomial eigenvectors, for any pair (P, Q) of polynomials, we 
have 

J P^Q)d^i = j Qh{P)dii. 

Since the coefficients gU and U are bounded on Ut with bounded coefficients, this 
identity may be extended to any pair (/, g) of smooth functions compactly supported 
in (not necessary with support in fi). Looking at this for smooth functions com¬ 


pactly supported in Ut leads to equation (2.3), which is equivalent to the symmetry 


property for such functions. Furthermore, applying this symmetry property to a pair 
of smooth function compactly supported in a neighborhood of a regular point of the 
boundary, and using Stokes formula, this implies in fact that, for any i — 1, • • • ,d, 
'^jg^^TXj = 0 at any point of the boundary, where (n^) is the normal vector to the 
boundary at that point. Therefore, this normal vector is in the kernel of the matrix 
g at any regular point of the boundary, which implies that det(^) = 0 at such a 
point. This gives item[^ 

We now know that the boundary is included in the algebraic set {det(^) = 0}, and 
we may look at the reduced equation for it, say Pi • • • = P = 0. Let x be a regular 

point of the boundary and V a neighborhood of it such that dUt H V = {Pq = 0} H V, 
for some g = 1, • • • , fc. In V, the normal vector (n^) to the boundary is parallel to 
diPq, so that we also have for all i, g^djPq = 0 on {Pq — 0} H V. But g^djPq 
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is a polynomial, which vanishes in V on the zeros of Pq. This implies (since Pq is 
complex irreducible) that 

( 2 . 10 ) = 

3 


where Li^q is a polynomial, the degree of which is less than since deg(^^- g^^djPq) < 
deg{Pq) + Ui. Then, equation (2.8) is just a rephrasing of ( 2.10[ ). This gives it em[4} 
If we now apply equation (2.8) and look at the value of P given by formula (|2.^, 


we see that, when the measure is gai,--,ak^ 




i,k^ 


and therefore is a polynomial with deg(6^) < a^. 

Therefore, for every n G N, the associated operator maps T-Lji into Pn- Moreover, 


the boundary equation (2.8) shows that for any pair of smooth functions compactly 
supported in M^, for the associted operator Lq 




j /Lai,-" — J ^Lq,^ ... (/)(i/iQ,^ ... , 

and this in particular applies for polynomials. Therefore, the operator Lai^...^ak is 
symmetric on the finite dimensional space Pn^ and this allows to construct a basis 
of eigenvectors for Lai^...^ak made of orthogonal polynomials. This gives item|^ 
The last item[^ that we shall not use in this note, is more technical, and relies on 
the fact that, looking at equation (2.3), any density measure p is such that di log p is 


a rational function, with singularities concentrated on {det(^) = 0}, and degree —a^. 
We refer to [3], where the arguments are developed, and which furthermore provides a 
complete description of the possible measures in the case where the reduced equation 
of dO. is not {det(^) = 0}. 

The proof of the reverse part of Proposition |2.1| is just a rephrasing of that of 
item[5l 


From Proposition 2.1, the important data are the set (open subset of 


bounded with piecewise smooth boundary given by an algebraic reduced equation 
Pi ''' Pk — 0), and the operator P, given by polynomial functions (^^•^), elliptic in 
12, satisfying the degree condition!^ and the boundary equation]^ 

To fix the ideas, we provide a few definitions 

Definition 2.5. 

1 . A polynomial domain Q C is a bounded open set in with boundary 
512 ineluded in some algebraie surfaee with redueed equation {P{x) — 0}^ and 
sueh that there exists some valuation {ai, • • • ,a^} and some elliptie eo-metrie 
P = {g'^^) on 12 with deg{g'^^) < ai + aj satisfying the boundary equation (2.10). 


2. A polynomial system (12, P) is given by a polynomial domain 12 together with 
the assoeiated eo-metrie P. 
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3. A polynomial model is a triple where (ri,r) is a polynomial system, 

and pi is probability measure on with smooth density p sueh that r(x% log p) = 
Si, with deg{Si) < a^. 

By definition, to each polynomial domain corresponds at least one polynomial 
system (there may indeed be many different co-metrics T associated with the same 
domain see m)- Moreover, we saw that to any polynomial system corresponds 
many polynomial models. 


Remark 2.6. The valuation is not unique. Beyond the trivial ehange (ai, • • • , a^) 
(cai,--- ^coq), the same polynomial model (or system) may eorrespond to various 
valuations. We shall make no effort to provide the lowest ones sinee in general a 
good ehoiee is provided by a simple look at the eo-metrie T. 


In [3], a complete description of all polynomial models is provided when the 
chosen valuation is the natural one (we give this description in Section 10 at the 
end of the paper for completeness). This description relies in an essential way on 
the fact that for the natural valuation, the problem is affine invariant, that is that a 
polynomial domain Tt is transformed into another one through affine transformations. 
This allows for an analysis of the boundary equation, and to the classification of 
algebraic curves in the plane for which the boundary equation (2.8) has a non trivial 
solution, through the analysis of the singular points of the curve and its dual. 

This affine invariance is lost for other valuations, since an affine transformation 
no longer maps the set Tin of polynomials with degree at most n into itself. This 
paves the way for the construction of new models. In what follows, we shall mainly 
concentrate on the construction of polynomial systems in dimension 2. These two 
dimensional models also provides new 3-d models through the use of 2-fold covers 
of our 2-d models. These two fold covers already appear in [3]. But even for the 2-d 
models which already appear there {Qu and of Section]^ some two-fold 

coverings appear as new. The reason is that in [3], only the models with natural 
valuation are considered. Here, even though the 2-d models may be considered with 
the usual valuation, this is no longer the case for their coverings. 


3 Constructing polynomial systems 


A generic way for the construction of polynomial models in dimension d is to consider 
some other symmetric diffusion operator L (often in higher dimension) and look for 
functions (Xi, • • • ,X^) such that, setting X = (Xi, • • • ,X^) 


L(X,) = B\X), r(X„X,) = G^^‘(X), 


where B^ and 


are some smooth functions. Then according to formula (2.6), 


L($(X)) = L($)(X), 


where 

ij i 
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This new operator L has as reversible measure fi which is the image of the reversible 
measure /i of L under X. This is often a good way to identify image measures, 
through equation ( |2.3| ). Then, L corresponds to a new triple (fi,r,/i), where is 
the image X(Q), T = (G^^) and ft is the image of /i. 

Definition 3.1. 

1. When we have sueh funetions (Xi, • • • ,X^) sueh that r(X^,Xj) = we 

say that (Xi, • • • , X^) form a elosed system for T. 

2. If moreover L(X^) = B'^[X), we say that we have a elosed system for L. 

It may happen that for some specific polynomial model (fi, T, /i) and some func¬ 
tions X = (Xi, • • • ,X^), X is a closed system for T, but not for L. 

Now, if L itself maps polynomials with degree n into polynomials with degree n 
(say with the usual valuation), if Xi is a polynomial with degree a^, and if B'^{X) 
and G^-^(X) are polynomials in X, then L provide a next polynomial model with 
valuation (ai, • • • , a^). 

It may also happen that this transformation x X — (Xi, • • • ,X^) is a diffeo- 
morphism, in which case we do not distinguish between those two models. If this 
diffeomorphism and its inverse are given through polynomial transformations, and 
if both are polynomial systems or models, we say that these systems or models are 
isomorphic. It is not always easy to see when a model is an image of another one, 
or when they are isomorphic. 

Apart of one specific case (example 7 of Section [^, all the models which appears 
in [3] may be constructed either from the Euclidean Laplace operator in acting 
on function invariant under the symmetries of a regular lattice (examples 1, 6 and 
11 in Section [To|), or from the spherical Laplace operator on the unit sphere C 
acting on functions which are invariant under some finite subgroup of 0(3) (all the 
other models of Section[To|). Here, we shall explore in a systematic way all the models 
that one may construct from the finite subgroups of 0(3). This construction may 
be also carried in higher dimension letting the spherical Laplace operator on 
act on polynomials in (that is on the restriction to of such polynomials). 

The spherical Laplace operator on may be described through its action on 
linear forms. If e is any vector in the Euclidean space M^, we look at the associated 
linear form e* : x e • x, and more precisely to its restriction to the unit sphere, as 
a function \-^ M. Then, for the Laplace operator L^ and its associated carre du 
champ r^, we have 

(3.11) L^(e*) = -{d - l)e*, r^(e*, /*) = e • / - e*/*. 

Therefore, looking at the canonical basis (e^) of M^, we see that any polynomial 
in the variables Xi (= e*) is transformed under L^ into a polynomial with the same 
(natural) degree. Moreover, the spherical Laplace operator commutes with all the 
elements of 0{d). Then, if we are given any subgroup of 0{d) and if we look at 
the set of polynomials invariants under the group action, L^ will preserve this set. 
If we may describe some polynomial basis for these invariant polynomials, then we 
expect to get in such a way a closed system, and therefore construct new polynomial 
models. 
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4 Invariant polynomials 


The theory of invariant polynomials has a long history going back to D. Hilbert, E. 
Noether, etc. It now plays an important role in coding theory and combinatorics (see 
[7]). In what follows, we provide a brief account which is useful for the understanding 
of our construction method, reducing to the case of finite subgroups of 0{n). We 
refer to mi for further details. 

Given any finite subgroup G of 0(n), any element g ^ G acts on the set of linear 
functionals (xi, • • • , We may consider its action on homogeneous polynomials in 
the variables (xi, • • • , Xn) and look for invariant polynomials, that are homogeneous 
polynomials which are invariant under the group action. If one denotes by dn the 
dimension of the vector space of invariant polynomials with homogeneous degree 
n, then Molien’s formula allows to compute the Hilbert sum F{G^t) = 
through 



(4.12) 


Moreover, the set of invariant polynomials may be represented as follows. First, 
there exist n algebraically independent polynomials (^i, • • • ,^ 71)5 called primary in¬ 
variants, and some other invariant polynomials ( 771 ,* •• ^gk) (fhe number of them 
may depend on the choice of the 0 ^), called secondary invariants, such that any 
invariant may be written as 


k 


^0(^1, • • • • • • , 0 n ), 


where Pi are polynomials (in the variables (^i,*** Moreover, each 77 ^ satis¬ 

fies some monic polynomial equation in the variables 9 — ( 0 j), that is satisfies an 
algebraic identity of the form 


dT + dT + • • • + Qi.piiS) — 0 , 


where Qi^k{^) are polynomials in the variables (^i, • • • , 9n). These algebraic relations 
are called syzygies. 

Furthermore, there are only primary generators if and only if the group G is 
generated by reflections, that is when G is a Coxeter group. 

In order to construct polynomial systems, we then consider finite subgroups 
of 0 ( 77 .), compute their invariants (primary and secondary when they exist), look 
at their restriction to the unit sphere (that is consider those polynomials modulo 
~ !)• They are no longer homogeneous, and, since always be 

considered as a primary invariant, we may reduce to 77 . — 1 primary invariants, plus 
some number of secondary invariants. We then let the spherical Laplace operator act 
on them. Since the spherical Laplace operator commutes with rotations, it preserves 
the set of invariant polynomials. Moreover, it maps polynomials with degree n into 
polynomials with the same degree. Constructing such polynomial systems amounts 
then to choose some family (Ci, * *' 5 Cp) of invariants, and look for Cj), expecting 
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that it may be written as • • • ,(p) (that is to provide a closed system for F). 

Then, the extra condition on the degrees will be automatically satisfied, where the 
valuation is defined through = deg((^^). The difficulty then is to produce such a 
closed system (Q) of algebraically independent polynomials. When such happens, 
we produce a polynomial system which is an image of the starting Laplace operator. 

In all the examples in dimension 3, one may always chose 2 primary invariant 
(^ 1 ,^ 2 ) to produce a closed system (this is no longer true in higher dimension, see 


Section 11). Moreover, when one adds one secondary invariant 77 , we always obtain a 


closed system with 3 variables (0i,02, v) — (Ci? C 25 Cs)- Now, it turns out that, if one 
forgets about the algebraic relations Q{0i^62^r]) = 0, and consider the polynomials 

a polynomial co-metric in dimension 3, it provides a poly¬ 
nomial model on a bounded domain in which has the surface Q(Ci 5 C 2 ?C 3 ) = 0 
(the syzygy) as a part of its boundary. Although the first construction with just 
the primary invariants (^ 1 ,^ 2 ) is not surprising (all our groups are sub-groups of 
Coxeter groups), the second property (construction of 3-dimensional models from 
2-dimensional ones through the syzygies) remains quite mysterious. 

Let us show this phenomenon in dimension 1, on the simpler form of the cyclic 
group Zji acting in M^, as a rotation in the complex plane with angle 27vjn. Writing 
z = x + iy^ we may choose as primary invariant X — and secondary invariant 

Y = ^{z'^). We now restrict them to the unit circle and let the spherical Laplace 
operator act on it. Using formulas (3.11), or the computations provided at the 
beginning of Section one sees that 

r(X,X) = 71 ^( 1 -X^), L{X) = -n^X, 


and therefore it provides a closed system for L which corresponds (up to the factor 
n^), to the classical Jacobi operator on (—1,1). Now, if we add the variable T, we 
get again a closed system for F, with co-metric 


/F(X,X) F(X,y)^ 
\T{X,Y) F(X,y)^ 


n 


1 - -XY \ 
-XY 


The determinant of this matrix is 1 — and indeed we have X^ -h = 1 in 

A - X^ -XY 
^ -XY 1 - 

is a metric on the unit ball 12 = {1 — X^ — > 0} C which corresponds to the 


our model (this is the syzygy relating X and Y). But the metric 


model 2 in the 2-d polynomial models of [3] in Section 10 The various probability 
measures for this model have the form Ca{l—X‘^ — Y‘^YdXdY^ whith a > —1. When 
a = (d — 3)/2, for some integer d > 2, this corresponds to the image of the Laplace 
operator on through the projection (xi, • • • ^Xd^i) G 1 -^ {xi — X^X2 — Y). 

This measure concentrates when a ^ — 1 to the uniform measure on the boundary 
S^. The case that we just described as the image of the Laplace operator on 
corresponds in this model to a limiting case when d ^ 1. 

This is this phenomenon that will remain valid in dimension 2 in the examples 
described below, although we will not provide such simple geometric interpretation 
for the various 3 dimensional models constructed from the syzygies. 
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5 Finite subgroups of 0(3) 

In our context, we shall restrict to finite subgroups of 0(3). We first describe 
them, following |4]. There are only five (families of) finite subgroups of 50(3), 
described by F. Klein, corresponding to the cyclic, dihedral, tetrahedral, octahedral 
and icosahedral respectively, denoted in what follows as Cn, respectively. 

The groups T, 0,X correspond to the elements of 0(3) preserving respectively the 
tetrahedron, the octahedron or its dual the cube, the icosahedron or its dual the 
dodecahedron. 

The finite subgroups of 0(3) are described in two ways. The first class is obtained 
adding the central symmetry J : x —x to one of the subgroups of 50(3). If G is 
such a group, we denote Gj this new group, with \Gj\ — 2|0|. 

The second class is obtained by those groups G of SO (3) which contain a sub¬ 
group Gi of index 2. A new group denoted 0i|0 is obtained as OiU{ J^, g G 0\0i}. 
This provides the groups T|0, Cn\X>n^ Cn\C2n^ where in the case of the cyclic 

and dihedral groups, the structure of invariants may depend on the fact that n is 
odd or even. A complete table of Molien’s formulas is provided in [4] together with 
the associated list of invariants (with however some error in the secondary invariant 
for the group X ). 

In the following sections, we shall describe the various invariants, and provide 
the polynomial models which they produce, both in dimension 2 with the primary 
invariants, and then in dimension 3 with the use of the secondary ones and their 
syzygies. 

Among the subgroups of 0(3), the following are Coxeter groups: {n even) 

and Vn\V 2 n {n odd) ; CnlVn, for all n ; T|0, Oj and Xj. They will yield the 
primary invariants and hence a closed system and a model. Among them, some 
were known: those obtained from C2|T^2 (coaxial parabolas), (the cuspidal 

cubic with secant), T\0 (the swallow tail), Oj (the cuspidal cubic with tangent). 
But Vn^j {n even) and X)n\V 2 n (^ odd) for n larger yield new models involving 
Tchebychev polynomials, not very surprising ; and Xj yields a nice model with an 
angle based on tt/S whose existence has been the initial motivation of this work. 

Since each other subgroup of 0(3) is a subgroup of one of these Coxeter groups, 
we obtain the higher dimensional models by adding the secondary invariant as aux¬ 
iliary variable. In most examples, if the equation of the boundary of the two dimen- 
sionnal model yield by the Coxeter group is P(X, K) = 0, then the equation of the 
boundary of the three dimensionnal models are either of the form — P(X, Y) = 0, 
of the form X{Z^ - P(X, K)) = 0 or of the form - XP(X, T) = 0 and the bound¬ 
ary of the three dimensional domain is either a bounded two leaves cover of the two 
dimensional domain, either the same but bounded also by a plane. The case of the 
groups Cn\C2n OT is Special in that we have more than one secondary invariant; 
each of them yield a different three dimensional system. 

From now on, the operator F will always be the carre du champ operator of the 
sphere 
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6 Cyclic and dihedral groups 

Let (x, y, z) be the standard coordinate system in On the unit circle {z = 0}nS^, 
we choose n equidistant points (e^,i = The group acts on them by 

circular permutations, which consist of elements of SO{3) with vertical axis and 
angle 27vjn. 

We first consider the complex function Z — x-\-iy^ with its conjugate Z — x — iy^ 
and observe that 


r(z, z) = r(z, z) = r(z, z) = 2 - zz, 

and 

r{z^ z) = 1 — r(2:, z) = —zZ, r{z^ z) = —zZ. 

With this in hand, we set 

Xn = + z^), Yn = Q{Z^) = _ z^). 

Z Zi 

The 3 variables {z, T^) are linked by the relation = (1 — 

Then, the table 


lT{z,z)) r(z,X„) T{z,Yn) \ 

r= r(x„,x„) r(x„,y„) 

V r(y,,y„)/ 

is given by 


(6.13) 


X — z^ —nzXn 

nH{l - z^r-^ - X^J 


-nzYn 
-n^XnYn 

n\il-z^r-^-Yl 


We may chose 6i — z, 62 = X„ as primary invariants; these are the invariants 
of the Coxeter group C„|X>n- Hence, consider (0i,02) = {z,Xn)- From table (6.13), 


we see that they form a closed system for T. Let Ti be the extracted matrix corre¬ 
sponding to the two first lines and columns from T. 


Ti = 




—716162 


Up to the factor the determinant of this matrix is P(0i, 62) = (1 — 6\)'^ — 62^ 
and according to n being odd or even, it has 1 or 2 irreducible factors. Then, it is 
quite immediate to see that 

ri( 6 'i,logP) = -2n0i,ri(02,logP) = - 271 ^ 02 . 


It satisfies therefore the boundary equation. When n is odd, the set = 

{P{01,02) > 0} C is bounded and Ti) provides a polynomial system. When 

n = 2p is even, P = P1P2, where Pi = (1 — Op — 02, P 2 = (1 — 0iY + 02- The area 
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Figure 1: and 



Figure 2: 


to illustrate the even case 


plane with Pi(0i, O2) > 0, ^ 2 (^ 1 , ^ 2 ) > 0 and 0i G (0,1) has P 1 P 2 = 0 as reduced 
boundary equation, and is again a polynomial system. In this model, we 

may chose deg(0i) = 1 and deg(02) = n, which comes from the sphere interpretation, 
but we may observe that we may as well chose deg(0i) = 1, deg(02) = n — 1. Observe 
that these domains correspond to the disk if n = 1 (model 2 in Section 10) and to 
the double parabola if n = 2 (model 4 in Section [^, but are new as soon as n > 3. 

We now add the variable Yn = rj in the figure, which is our secondary invariant 
(observe that the roles of and are similar, and we may as well exchange them). 
This reflects the symmetries of the cyclic group We now have the co-metric 


r2 


/I - 0l -716162 

I n^((l — 6l)^~^ 


—nOiTj 

62) —n?62ri 

n^((l — 6\)'^~^ 


whose determinant factorizes as n^(l — — 

last factor 


p(ei,e2,v) = {i-olr-el-v^ 




rf). Note that the 


reflects the syzygy relating r] to (0i,02)- It is not a surprise that this determinant 
vanishes identically, since F is the Gramm matrix of the gradients (on the sphere) of 
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three functions, and the range of these 3 gradients is at most 2. Observe also that 
this syzygy, which will appear in the boundary of the 3-d system, may be written as 
^2) ~ where P is the boundary equation of the corresponding 2-d system. 
But now consider r2 as a co-metric in on the bounded domain = {|0i| < 
1, P(0i, ^2, '^) > 0} C M^, which has indeed again reduced equation P(0i, ^2, '^) = 0. 
We may check that 

r2(0i,log(P)) = -2n0i,r2(02,log(P)) = -2n202,r2(r?,log(P)) = -2n\ 

so that indeed (i72”\r2) is a polynomial system in with degrees deg(0i) = 
l,deg(02) = n,deg(r?) = n. 



Figure 3: Here, the domains 



and 


We now study the case where the groups contain the central symmetry. This 
corresponds to the new system of primary invariants {6\ = z'^, 62 = W) associated 
with the Coxeter group T>n,J (n even) or Dn\'D2n (n odd). We get 

fidiil-di) -2n0i02 A 

which, up to a constant, has determinant P(0i, ^2) = 6^2) ~ 

It has three irreducible components when n is even and 2 when n is odd. Once again 

r3(0i,log(0i)) =4(1 -0i),r3(02,log(0i)) = -2n02, 


and, 

r3(6'i,log(Fi)) = -4n6'i,r3(6'2,log(Pi)) = - 271 ^ 02 , 

so that the domain ^3 = {0\ G (0, 1),P\{0\,02) > 0}, which has reduced bound¬ 
ary equation 0iPi{0ii02) = 0 is such that (H3,r3) is a polynomial system. 
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Figure 4: The domains and (for Cs^j and C^^j) 



Figure 5: The domain 




(4) 

3 


Let us now add the secondary invariant rf = Yn^ to treat the group Cn^j when n 
is even and Cn\C2n when n is odd. We get a new co-metric 


r4 = 


4(9i(l-6>i) -2n6>i6>2 

n^((l — 


—2n6ir] 

—ri^02rj 

n2((l - - rj^)^ 


The determinant of this matrix factorizes as 


n\i-err-Hi{{i-err-el-rf). 

The factor ^2, '^) = (1 ~ ^1)^ ~^2~d^ represents the relation between 77, 0i, 62 

(the syzygy). The two factors Fi(0i, ^2,'^) = (1 ~ ~ ^2 ~ satisfy the 

boundary equations 


r4(0i,log(0i)) = 4(1 -0i),r4(02,log(0i)) = -2n02,r4(r?,log(0i)) = -2nr] 
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r 4 ( 6 'i,log(Pi)) = - 4 n 6 'i,r 4 ( 6 ' 2 ,log(Pi)) = - 2 n^ 6 ' 2 , r 4 ( 77 , log(Pi)) = -2n^r?. 

(This is not true for the factor 1 — 0i). The domain f 24 C defined by 0i G 
( 0 , 1 ), Pi( 0 i, 02 ,'^) > 0 , which has reduced boundary equation 01 ^ 1 ( 01 , 02 , 0 ) = 0 , is 
therefore such that (fi 4 , r 4 ) is a polynomial system. Observe that the syzygy, which 
appears in one component of the boundary of fi 4 , may be written as P( 0 i, 02 ) — 0 ^, 
where P is one of the components of the boundary of the corresponding 2-d domain 



Figure 6 : The domain Q 


(3) 

4 


We may now consider the dihedral group which amounts to add to the sym¬ 
metries of Cn the transformation (x, z) (x, —y^ —z). It has as primary invariants 
01 = z ^,02 = Xn as before (corresponding to the 2 -d polynomial system 
but now the secondary invariant is 0 = zY^. The new co-metric in dimension 3 is 
then 

/40i(l-0i) -2n0i02 -27/((n-h 1)01 - 1) 

Fs = n 2 ((l - 0 i)^-i - 02 ) -n{n + 1)020 

\ (1 — 0 i)’^“^(l -h (^2 — l) 0 i) — 02 — (n -h 1 )^ 0 ^ 

The determinant of this metric factorizes as 

4n2(0i(l - OiT - Oiel - r? 2 )((l - 0 i)--i((n 2 - l) 0 i - 1) - Ol) = Ar?PiP2, 
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where Pi(Oi, ^ 2 , d) = ^i(l ~ ~ ^ 1^2 ~ syzygy which relates 77 to ( 0 i, 02 )- 

Observe once again the relation between this syzygy and the boundary equation of 
the corresponding 2-d domain fis. 

Once again, we have 


r5(0i,log(Pi)) = 4(l-(n+l)0i),r5(02,log(Pl)) = -2n(n+l)02,r5O,log(Pl)) = -2{n+lfv 


while the boundary equation is not satisfied for P2. In the domain fls delimited 


by 0i G (0,1), Pi > 0 is a bounded domain with reduced boundary equation Pi = 0, 
and (f] 5 ,r 5 ) provides a 3-dimensional polynomial system. 



Figure 7: The domain O 5 for n = 3 


The groups T)n\P2n for even or for n even have primary invariants {z‘^ = 
01 , = 02) and secondary invariant 771 = zY^ or 772 = zXn- However the groups 

Cn\C2n even) or Cn^j {n odd) have the same primary invariants and as secondary 
invariants zYn^ X^Y^. It may be worth to observe that (z, X^, XnYn) is another 
form of the invariants for C2n, since X 2 n = X^ — Y^ — 2X^ — (1 — z^)^, and 
Y2n = 2 X^ 17 ,. and therefore they do not provide any new model (although they 
provide them under another form). 

We first choose 0i = 2 ;^, 02 = X^ (for which we already know that it corresponds 
to {z‘^^X2n) through a change of variables. We then get a co-metric 



0l(l — 0i) —710102 

n202((l-0i)^-i-02) 
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which corresponds to a 2-d domain fig with boundary reduced equation 01^2((1 ~ 
9i)^ — 62 ) = 0, which is isomorphic to the domain when changing n into 2 n (this 
model has 3 irreducible components in its boundary equation, and two of them may 
be reduced to a line, providing then a simpler form). 



Figure 8: The domain VLq for n = 3 


We may first add the secondary invariant rji = zY . We get a co-metric 

/40i(l - Oi) - 4 n 0 i 02 2771(1 -{n + l)0i) 

r7= 4n202((l- 02 ) -2n(n + 1)0201 

\ (1 - 0 i)"'“^(l + {v? — l) 0 i) - (n + 1)^01 - 02 

The syzygy relation between (^ 1 ,^ 25 '^) niay be written as 

P(01,02,0i) = 0i((1-0i)--02)-02, 

once again of the form Q(0i, ^ 2 ) — where Q appears in the boundary equation of 
the corresponding 2 -d domain fig. 

One may check that for this co-metric, O2P divides det(M), and moreover that 
r 7 ( 0 i,log(P)) = 4(l-(n+l)0i), r 7 ( 02 ,log(P)) = - 4 n(n+l) 02 , r 7 ( 0 i,log(P)) = - 2 (n+l)^ 0 i. 
and also 

r 7 ( 0 i,log( 02 )) = -4n0i, r 7 ( 02 ,log( 02 )) = 4 n^((l- 0 i)’"“^- 02 ,r 7 ( 0 ,log( 02 )) = - 2 n(n+l) 0 i. 

This provides a 3-d domain C with boundary reduced equation 02^(01,02, '^ 1 ) = 

0 , such that (fl 7 ,r 7 ) is again a polynomial system. 

Adding the variable 772 = zX^ instead of 771 = zT, to = z^, 02 = leads to 
the co-metric 
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Figure 9: The domain for n = 3 


/40i(l - Oi) -4n0i02 2772(1 - (n + l)0i) \ 

Fs = I 4 n^ 02 ((l - - 62) 2777/2(^(1 - 9 iY~^ - (77 + 1 ) 02 ) I 

\ 77^0l(l - 9iY~^ +92- {n+ 1)^02/ 

The determinant of this matrix has 3 factors, 2 of them being Pi — 9i92 — rj^ and 
P 2 — 92 — {l — 9i)'^. Pi(0i, 02, 02 ) is the syzygy relating 7/2 to (0i, 02). It is still of the 
form Q (01, 02) —'^2’ where Q appears in the boundary equation of the corresponding 
2-d domain fig. 

Now, once again, we have 

rr8(0i,log(Pi))= 4(1 -(77+1)01), 

\ r8(02,log(Pi)) = 4n((l - - (n + 1)02), 

[r8(??2, log(Pi)) = -2(n + 1)^772, 


'r8(0i,log(P2)) = -4n0i, 

^ r8(^2, log(P2)) = —4n^02, 
Jsim, log(P 2 )) = -2n(n + l)r/ 2 - 
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The third factor of the determinant does not satisfy the boundary equation. The 
domain ftg with boundary reduced equation P 1 P 2 = 0 provides a polynomial system 

(f^8,r8). 



Figure 10: The domain Qg for n = 3 


We now may add instead the secondary invariant 773 = X^Yn- As already ob¬ 
served, the new system (z^ = ^ijX^ = — ^nYn) is isomorphic to the system 

{9i = = X2n^d — Y2n) described by the metric r 4 and the domain Ob¬ 

serve however that in this presentation, the reduced boundary of the domain is 
0i(02(l — Oi)'^ ~ ^2 ~ Vs) — O 5 i^he second factor being the syzygy. 

Finally, one may check that adding two of the secondary invariants rji to 0i — ^ 

62 — X^ will not provide any closed system for F. The system (^i, ^ 2 , Vi^V 2 ^ Vs) pro¬ 
vides a closed 5 dimensional F operator, but the determinant of the metric vanishes 
(in M^) and there does not seem to be any polynomial system associated with it. 

7 Tetrahedron and cube/octahedron 

It makes sense, as we shall see, to treat jointly these cases. The groups T, O corre¬ 
spond to the elements of SO(3) preserving respectively the tetrahedron, the octahe¬ 
dron or its dual the cube. They have respective cardinality 12 and 24. Adding the 
central symmetry J : x —x to each of them we obtain 7j and Oj. Observe that 
the first one does not preserve the tetrahedron, while the second one does preserves 
the cube. We also consider the group T\0 which can be obtained by adding a plane 
symmetry with respect to the plane symmetry axes of the tetrahedron and which 
preserves the tetrahedron regardless of orientation. 
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Figure 11: The domain Og for n = 3. 


They are related by the following inclusions diagram : 

r 

Tj O T\0 


Oj 


Let {x,y,z) be the standard coordinate system in We put the cube centered 
at the origin and with faces parallel to the coordinate planes. We put the tetrahedron 
with edges on the diagonal of the cube. We consider the polynomials 

O 3 = xyz ,04 = x^ + y^ + z^,Oq = - z^){z‘^ - 

which will play the same role as the one played by (z^X^Yn) in the previous one as 
basic blocks to construct all the invariants for the various groups concerned in this 
section. 

We first compute : 


/r(03,03) r(03,04) r(03,06)\ 
r= r(04,04) r(04,06) 

V r(06,06)/ 


which is 


/(I - 04)72 - 90| 403(1 - 3O4) 

8(60i + 3 O 4 - 1 - 20|) 


—I 8 O 3 O 6 
806 ( 2 - 304 ) 

-540|04 + 180| - 30| + 4 O 4 - 1 - 360| 
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We consider the primary invariants of the Coxeter group T\0 given by {9 1 ^ 62 ) = 
( 03 , 04 ). The determinant of the submatrix Tn given by the first two rows and 
columns is 


P{01,O2) = -lO80f + 200^ + 261 - 5^1 + 402 - 360^02- 

It provides a domain fin with boundary P{X^Y) — 0. This corresponds to the 
model of the swallow tail (example 10 in Section [Tol). 



Figure 12: The domain flu : swallow tail. 


For the group 7^, we can choose 77 = Oe as secondary invariant. It is alge¬ 
braically related to (^ 1 ,^ 2 ) = ( 03 , 04 ) through 

= PiOl, 62 ). 


If we write the matrix 


/r(0i,0i) r(0i,02) r(0i,r?)\ 

ri2= r(02,02) r(02,r?) , 

\ r(0,0)/ 


we get 

/ _90i2_ 02 / 2 + 1/2 -120102 + 401 

48 0^- 1602^ + 24 02 -8 




— 18 Biri 
— 24927 ] + 1677 

540^02 + 180 ^ - 302 ^ - 367/2 + 4(92 
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The determinant of this matrix factorizes in 

= 4 (3 02 - 1) (18 01^ + 02 - 1) (108 01^ + 36 01^02 - 2 02^ - 20 01^ + 5 02^ + 47/2 - 402 + 1) 
It turns out that the factor 

-P3(01,02,0) = lO80f+ 36 0^02 - 20^- 20 0? +50| +47/2- 402 + 1 

satisfies 

rri2(0i,logP3) = —3601, 

Ti2(02,log(P3)) =-4802 + 32, 

[ri2(7/,l0g(P3)) = -727/, 

SO that this provides a new polynomial model (rii2,ri2) in dimension 3 . (The 
boundary equation is not satisfied for the two other factors.) 

We can check that the complementary of the surface Ps{X^Y^ Z) = 0 has one 
bounded component in and that the determinant does not vanish inside this 
component. 



Figure 13: Surface Z) = 0, bounding Yti 2 . 


We observe that Oj is also a Coxeter group. We can take as primary invariants 
( 01 , 02 ) = ( 0 ^, 04 ) for Oj. We get 

^ /r(0i,0i) =40i((1 - 02)/2 - 90i) r(0i, 02 ) = 80i(1 - 302) \ 

V r(02,02) = 16(301 + |02-1/2-02)^ 
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whose determinant is given by, 

02) = 160i(-lO80^ + 2001 + 20l - 50i + 402 - 1 - 360102) 


Observe that Q(X^,Y) = X^P(X,Y). We recognize that the 
domain is the cuspidal cubic with tangent (model 9 in Section 


boundary of the 

To). 



Figure 14: The domain O 13 : cuspidal cubic with tangent. 

For the group 7j, we may add the secondary invariant ig — Oq. We obtain the 
co-metric ri 4 : 

■ -36 01^ - 20102 + 201 -24 0i02 + 8 0i -360ir? 

4801-1602^ + 2402 - 8 -24027?+ 16r/ 

-54 0102 + 1801 -3 02^-367?2 +402 -1 

The determinant of this matrix factorizes as 

160i( 302 - l)(180i + 02 - l)(-20i + 1080? + 360102 + 50? + 47?^ - 2O0i - 402 + 1). 

Only the two factors 0i and 

^3(01) 02) v) — (—20? + 1080? + 360102 + 50? + 4??^ — 2O0i — 402 + 1) 
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satisfy the boundary equation, with 

'ri4(0i,log(0i)) = -3602-202 + 2, 

< ri4(02,log(0i)) = —2402 + 8, 

^ri4(r7,log(0i)) = -36r], 

and 

'ri4(0i,log(Q3)) = -720i, 

< ri4(02,log(Q3)) =-4802 + 32, 

,ri4(03,log((33)) = —7303. 

We observe that this factor writes r]^ — Q(0i,02). It appears that the two com¬ 
ponents bound a domain I7i4 on which the other factors do not vanish. 



Figure 15: Surface XQ 3 (X,Y, Z) = 0 bounding the domain Vti 4 ^ 


Finally, for the group O, we use (0i,02,??) = (0^,04 , OsOq). We compute 
the Fis matrix 

■-36 01^- 20102 + 201 -24 0102 + 8 01 -540ir? - 02/? + r? ' 

48 01 - 16 02^ + 24 02 - 8 -36r/02 + 20r] 

^ 3,3 
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with 

(^3,3 — 8177 ^- ^ —^+9/2301^2+3/201+3/2^- +5/802~l/3* 

2 8 8 

The determinant of the matrix factorizes as 2QiQ2^ with 

Qi^ 2 02^ + 324 01^02 - 12 0102^ - 7 02^ - 36 0i^ + 24 0i02 + 9 02^ - 12 0i - 5 02 + 1 

Q 2 = 2 0102^ + 108 01^ + 36 01^02 + 5 0102^ - 20 0i^ - 4 0i02 + 4 77 ^ + 0i. 

Only Q 2 satisfies the boundary equation, with 

rri5(0i,log(Q2)) = 2-lO80i-202 
lri5(02,log(Q2)) = 40 - 7202 
[ri5(r/,log(Q2)) = -162r/ 

We observe that (OsOe)^ = 03 ^( 03 , 04 ) so that rj^ — 01^(01)02) =: -R(0i)02)- We 
get then a new domain 14 15 with boundary 

Q 2 {Oi, 02 , rf) := rf - i?( 0 i, 02 ) = 0 . 

We may observe that in all these cases, the 3-dimensional domains can be cut by the 
plane {Z — 0} to get a 2-dimensional domain. In the first case, we get the swallow 
tail, in the other two cases, we get the cuspidal cubic with tangent. They provide 
various two fold coverings of these two dimensional models. It could be interesting 
to investigate the shape of the singularities of the boundaries of these 3-dimensional 
domains. 

8 Dodecahedron / Icosahedron 

We finish by the study of the groups X and Xj of cardinality 60 and 120. The 
computations are a bit more painful but the idea is always the same. Let us introduce 
as before three new building blocks. With c = (1 + \/5)/2, 

/e = {(?x^ — y^){c^y^ — z^)(c^z^ — x^), 

Iio = (^ + y + z)(-x + y + z)(x - y + z)(x + y - z) 

(c~^X^ — — C^Z^))(c~^Z^ — C^X^), 

7 i5 = xyz{cx + c~^y + z)(—cx + c~^y + z)((cx — c~^y + z){cx + c~^y — z) 

{x + cy + c~^z){—x + cy + c~^z){{x — cy + c~^z){x + cy — c~^z) 

{c~^x + y + cz){—c~^x + y + cz){c~^x — y + cz){c~^x + y — cz). 

We will use the primary invariants ( 0 i, 02 ) = (/6,7io) and the secondary invariant 
77 = / 15 . Let us do the computations directly with all the invariants in this family. 

With ( 01 , 02 ,'/?) = (76 ,/105 7 i 5 ), we get for the co-metric 

/r( 0 i, 0 i) r( 0 i, 02 ) r(0i,77)\ 

r= r( 02 , 02 ) r(02,77) 

V r(r?,r/)/ 
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Figure 16: Surface Z) = 0 bounding the domain 


with 

'r{0i,ei) = -36er - (V5 + 2) (701 + 502 + 2^502), 
r(0i, 02) = (40 - 16 Vs) 0V + (3 Vs + 6) 02 + 01 Vs - 60 01 02, 
r(02,02) = (7296 - 3264 Vs) 0V + (96 VS - 240) 0i 02 

+ (-432 + 192 Vs) 0V - 5 V502 + (6 - 3 Vs) 0i - 100 02^ 

< r( 0 i,V = -90017 /-2 (Vs+ 2 ) 7? 
r(02,7?) = -150 02 7?+ 7? ((- 100 + 40 Vs) 01 -2 Vs) 
r(7/,7/) = -225 7/2 - 1/4 (-161 + 72 Vs) 

(13 01 Vs + 45 V502 + 45 01^ + 4 Vs + 26 01 + 100 02 + 9) 

(30^1^2 Vs-3 01 Vs-9 V502 - 19 0V+ 75 01 02 - 6 01 - 20 02) 

For the Coxeter group Ij we restrict our attention to (0i, 02)- Up to some factor, 
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the determinant of the sub-matrix T 21 corresponding to (^ 1 ,^ 2 ) is S{ 0 i^ 62 )^ with 

5(01, 02) = 688 VEef + 6480 \/50?02 + 1728 Of + 364 efVb + 6042 \/50?02 + 23400 VEOi 0l 

+17050 \fl0l + 1376 0\ + 14400 0\02 + 68 0? + 1288 0i 02 + 1220 ^50^ 

-19520 + 819 0\ + 13515 0?02 + 52325 0i 0\ + 38125 0\ + 152 0\ + 2880 0i 02 

+2728 0i 

This is a new 2-dimensional domain ^ 21 - 



Figure 17: Domain D 21 bounded by S{X^Y) — 0. 

For the direct subgroup X, the determinant of r 22 factorizes (up to some 
constant) into 818283 ^ with 

5i(0i,02,0) = 5'(0i,02)-43648r?2, 

32 ( 01 , 02 ) = 13 01 \/5 + 45 V 5 O 2 + 45 0? + 4 \/5 + 26 01 + 100 O 2 + 9, 
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and, 


-53(01,02) = 30 01 02 V5 - 3 01 V5 - 9 ^502 - 19 0? + 75 0i 02 - 6 0i - 20 02- 

^ 2 ,'^) is the syzygy relating r] to 9i and 02 , and the polynomial iSi satisfies 
the boundary condition for the co-metric r 22 , which once again provides a new 
polynomial system with domain Q 22 in dimension 3, since the boundary conditions 
are satisfied: 

'r22(0i,log(5i)) = -4^5 - 8 - 18001, 

< r22(02,log(5i)) = 4(2^5-5)((30^+ 75)02 + 1001 + 2 + ^5), 
V22(r?,log(5i)) = -450r?. 

Observe that we may as well rescale rj to get a simpler domain tf = S'(0i, 02). 



Figure 18: Surface Si{X,Y, Z) = 0 bounding the domain 022- 
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9 Summary 

We summarize here the models detailed above. For the dihedral family, we let 

Hn(X,Y) = (l-Xr-Y. 


Group 

01 

O 2 

V 

n 

Boundary 

Picture 

C W 

Z 

Xn 



Hn{X\Y^)=0 

0 /V 

Cn 

Z 

Xn 

Yn 

S!<”' 

Hn{X‘^,Y^) - = 0 

^ 4 

'Y’n,J: 'Un\D2n 


Xn 


S!<“) 

XHniX,Y^) = 0 


\ 

\ 

/ 

/ 

CnJ 5 C 77 , |C 2 n 


Xn 

Yn 


X{Hn{X,Y^)-Z^)^0 


Vn 

z^ 

Xn 

zYn 

Sif 

XHn{X,Y‘^) - Z^ ^0 


Y>2n,J 





XYHn{X,Y)^0 




'Un\U2n, 'Un,J 

Cn |f' 2 n) Cfi^J 



zY 

^-L n 


XHn{X,Y)- Z^ ^0 


Cn 1^2717 Cn^J 

^2 

y‘^ 

zXn 

nP 

{XY-Z‘^)Hn{X,Y)^0 


Cn 1^277,7 Cn^J 

^2 

y‘^ 

XnYn 


X{YHn{X,Y)-Z‘^)^0 

y 


When there are two groups, the first one is for n odd, the second is for n even. 
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For the tetrahedron / cube / octahedron family, we let H(X,Y) = 108X^ — 
20X — 2Y^ + bY^ — 4:Y + 36Xy, and for the dodecahedron / icosahedron family, we 
set S be as defined in Section [H 


G 

0i 

O 2 

V 

n 

Boundary 

Picture 

r\o 

03 

O 4 


f2ii 

H{X^,Y) =0 

A 

T 

03 

O 4 

Oe 


-4Z2 = 0 


Oj 

Ol 

O 4 


f^l3 

XH{X,Y) = 0 


Tj 

01 

O 4 

Oe 

f2l4 

X{H{X,Y) = 0 


O 

Ol 

O 4 

OsOq 

i7i5 

Z‘^ - XH{X,Y) = 0 


Tj 

Oe 

Oio 



5(X,F) = 0 


X 

Oe 

Oio 

Oi5 

^22 

= S{X, Y) 

> 


Covers If is a subgroup of G, then polynomials that are invariant by G are also 
invariant by H. It follows that we can express the G invariants (primary and sec¬ 
ondary) as polynomials in terms of the primary and secondary i7-invariants (modulo 
+ y‘^ + — 1). These polynomials define a mapping from the i7-domain onto the 

G-domain which is a /i-covering, h being the index of H in G. 

Let us see what happens for a few examples. This is specially easy for the 2- 
covers. We go from 3D to 2D models by simple projection (forgetting variable z). 
In other cases the effect is that of adding a new symmetry for instance : 

{x,y) iH- {x^,y) 
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(x,y,z) iH- (x^,y,xz) 

{x,y,z) 1 -^ {x‘^,y,z) 

The case of C 3 IP 3 as subgroup of T|0 is a little bit more tedious because we did 
not choose the same coordinates for the representations. We can see what happens 
taking for C 3 IP 3 the variables {x + y + z)/^/3 and xyz to have a common invariant. 
It is a linear computation to get a, 6 , c and d such that 

x^+ 7 /^+ 2 ;^ = a{x-\-y-\-z)^+b{x-\-y-\-z)‘^{x‘^-\-y‘^+z‘^)+c{x-\-y-\-z)xyz-\-d{x‘^-\-y‘^-\-z‘^)‘^. 
We obtain a map of the form : 

fin 

{x^y) {y ^ ax"^-\-bx^^-\-cxy-\-d). 

10 The bounded two dimensional models of [3] 


We provide here for completeness the complete list of models in dimension 2 de¬ 
scribed in j3]. With the restriction that the valuation is the usual one, they are the 
only ones which may occur up to affine transformations. We indicate the (scalar) 
curvature when it is constant (+ when it is a positive constant, 0 otherwise). When 
no curvature is indicated, it comes from the fact that the metric is not unique (mod¬ 
els 2 and 3), in which case there exist at least one metric for which the curvature 
is constant and positive), or it is not constant (model 7). Up to isomorphism, one 
may replace one parabola by an horizontal line in model 4, so that this changes the 
degree in the boundary (in this particular case however, the co-metric is no longer 
unique). 

Up to isomorphism, we have 

f (2) - (3) fi^^\ (4) ~ fif^ ~ 

1 (5) ~ fi^^^ ~ fi^‘^\ ( 8 ) ~ fi^^\ (9) ~ fii3, ( 10 ) ~ fill 
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tt 

Curv. 

d(Q) 

Boundary 

Picture 

1 

0 

4 

(i-x2)(i-r2) ^0 

□ 


2 


2 

1 - x2 - y2 ^ 0 



3 


3 

XY{l-X-Y) = 0 

A 

4 

+ 

4,3 

(1 -X2)2 -y2 ^0 

0 

5 

+ 

4 

Y{1- X){X‘^ -Y) 

/ 

7 

6 

0 

4 

{Y - X2)((y + 1)2 - 4X2) ^ Q 

A 

7 


3 

y2 _x2(i -X) = 0 

c>< 

8 

+ 

4 

(y2 _x3)(x - 1 ) = 0 



9 

+ 

4 

(y2 - x3)(2(y -1) - 3(x -1)) = 0 

/ 

10 

+ 

4 

4X2 _ 27X4 + ley - I28y2 - 144X2y + 256y3 = 0 


11 

0 

4 

(x 2 + y 2)2 + i 8 (x 2 + y 2 ) - 8x3 ^ 24xy2 - 27 = o 

A 


11 Further remarks 

In the various models presented here, it happens that the primary invariants provide 
a closed system. The reason is that all these groups are subgroups of finite Coxeter 
groups, for which the invariants are our primary invariants. It is not true that this 
is always the case. Here is an example provided by Y. Cornulier of a group in 
dimension 4 which is not a subgroup of a finite Coxeter group. Let Mp the matrix 
of a rotation with angle 27t/p in M^, where p is an odd prime. Let I 2 be the 2x2 
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identity matrix. Then, we consider the group generated by 










This group has elements, A^ 2 ^, 0 < ki, k 2 < p 

Mollien’s formula (4.12), the Hilbert sum is easily computed 


1 , and, thanks to 


m- 


i+tp + 2tp+^ + 2fp+ep+^ + t^p+^ 

(1 - f2)(l - t^)(l - t^P)[l - tP) 


This leads to the description of primary and secondary invariants when restricted 
on the unit sphere in Following Section]^ we identify C^, and for a pair 

{zi,Z 2 ), consider zj = Xj + iYj and Rj = \zj\‘^, j = 1,2. Then we chose 


9 i = Xi + X 2,92 = X1X2,93 = R1R2 


as primary invariants, and the secondary invariants may be chosen as 

Vi^Yi + Y2,m = {Yi - Y2 ){Ri - R 2 ),m = (^1 - X2){Ri - R2), 

r?4 = T 1 T 2 , V5 = {Xi - X 2 )iYi - Y2),m = (^lU - X 2 Y 2 ){Ri - R2),V7 = V3V4 

It turns out that 9\,92,93 is not closed for F, where F is the square field operator 
on the unit sphere in R^. For example 


r(0i, 02) =y(Tp_i-20102), 

where = R 1 X 2 + R^Xi, and Tp-i may be expressed as 

Tp-i — 0i(3i(03) +Ti ( 32 (^ 3 ), 

where Qi are polynomials, and Ti = [Oi + rj^)/2^ so that r(0i, 02 ) is not a polynomial 
of ( 01 , 02 , 03)- 

In this example, one may observe that indeed (0i, 02, 03 , '^ 3 ) form a closed system 
for r, but this does not provide any model in (the boundary equation is not 
satisfied). 

A final remark is that our various 3-d models constructed from 2-d one could 
appear as provided by Coxeter groups in dimension 4. If such would be the case, 
the natural Ricci curvature carried by the associated cometric T would be constant 
(since it would locally be the spherical co-metric seen through a diffeomorphism). 
One may easily check that this is not the case. 
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